Classical memory effects on spin dynamics in two-dimensional systems 
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We discuss classical dynamics of electron spin in two-dimensional semiconductors with a spin-split 
spectrum. We focus on a special case, when spin-orbit induced random magnetic field is directed 
along a fixed axis. This case is realized in III-V-based quantum wells grown in [110] direction and 
also in [100]-grown quantum wells with equal strength of Dresselhaus and Bychkov-Rashba spin-orbit 
couplings. We show that in such wells the long-time spin dynamics is determined by non-Markovian 
memory effects. Due to these effects the non-exponential tail 1/t 2 appears in the spin polarization. 
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PACS numbers: 71.70Ej, 72.25.Dc, 73.23.-b, 73.63.-b 

Continuous reduction of the device sizes in last decades 
has initiated active research of transport, optical and 
spin-dependent properties of low-dimensional nanostruc- 
tures. In recent years, it was clearly understood that not 
only quantum but also purely classical phenomena might 
lead to rich physics in such structures. In particular, 
a number of non-trivial transport phenomena, such as 
magnetic- field- induced classical localization 0,12], high- 
field negative 0, 0, @, 0, and positive [|| maffnetoresis- 
tance, low- field anomalous magnetoresistance j3,@i zero- 
frequency conductivity anomal y fej , and non-Lorentzian 
shape of cyclotron resonance [lfj might be realized in 
two-dimensional (2D) disordered systems. All these phe- 
nomena arise due to classical non-Markovian memory ef- 
fects which are neglected in the Drude-Boltzmann ap- 
proach. The strength of these effects is governed by a 
classical parameter d/l (d is the characteristic scale of 
the disorder and / is the transport scattering length). 
Since the role of quantum effects is characterized by a 
parameter X/l (A is the electron wavelength), the clas- 
sical effects might dominate in systems with long-range 
disorder, where d^$> X. 

Usually, classical memory effects slow down relaxation 
processes leading to non-exponential decay of correla- 
tion functions. In particular, the velocity autocorrela- 
tion function in a 2D disordered system has a power tail 

dm 



Mt)v(o)) 



= -c 



f 2 



t > T, 



(1) 



in contrast to exponential decay exp(— t/r) predicted by 
the Boltzmann equation. Here r = 1/vf is the trans- 
port scattering time, vf is the Fermi velocity and C is 
the coefficient which depends on the type of disorder: 
C = 2d/37rZ for the Lorentz gas model, where electrons 
scatter on hard disks of radius d randomly distributed 
in a 2D plane with concentration n (nd 2 < 1) and 
C ~ d 2 /l 2 for the scattering on the smooth random 
potential with a characteristic scale d. Physically, this 
long-lived tail is due to "non-Markovian memory" spe- 
cific for diffusive returns to the same scattering center 



[lij ] (see also recent discussion in Refs. [Uli^). 

In spite of the large number of publications devoted to 
the study of non-Markovian transport phenomena, the 
role of memory effects in spin dynamics is not well un- 
derstood. In this paper we discuss the slow down of the 
spin relaxation in 2D systems due to the non-Markovian 
memory. This effect is of particular interest for new 
rapidly growing branch of semiconductor physics, spin- 
tronics. The main goal of spintronics is the development 
of novel electronic devices that exploit the electron charge 
and spin on equal footing fl4j . For effective functioning 
of such devices, the lifetime of the non-equilibrium spin 
must be long compared to the device operation time. In 
III-V-based semiconductor nanostructures, this require- 
ment is not easy to satisfy, since in such structures the 
spin polarization relaxes rapidly due to Dyakonov-Perel 
(DP) spin relaxation mechanism [151] . This mechanism 
predicts the exponential relaxation of non-equilibrium 
spin with a certain characteristic time ts- At small tem- 
peratures, this relaxation m ight slow down due to quan- 
tum interference effects However, with increas- 
ing temperature, the interference effects are suppressed 
by inelastic scattering. Here we show that the classi- 
cal non-Markovian effects which are not very sensitive 
to the temperature might lead to long-lived 1/t 2 non- 
exponential tail in the spin polarization in analogy with 
velocity relaxation described by Eq. (JJJ . 

The DP mechanism is based on the classical picture of 
the angular spin diffusion in random magnetic field in- 
duced by spin-orbit coupling. In 2D systems, the corre- 
sponding spin-relaxation time t$ is inversely proportional 
to the momentum relaxation time r: l/rg ~ (Wp)r |l8| . 
Here w p is the frequency of spin precession in a random 
magnetic field, p is the electron momentum, and angu- 
lar brackets denote averaging over momentum directions 
(for p w Pf)- As a consequence, in high-mobility struc- 
tures which are most promising for device applications, 
ts is especially short. However, in some special cases, the 
relaxation of one of the spin components can be rather 
slow even in a system with hi gh m obility. In particular, 
a number of recent researches ]H EoL Ell El I2I El] are 
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devoted to GaAs symmetric quantum wells (QW) grown 
in [110] direction. In such structures, to p is perpendicular 
to the QW plane 0] and depends on one component of 
the in-plane momentum (say x-component) 



ap x z. 



(2) 



Here z is the unit vector normal to the well plane and 
a characterizes the strength of the spin-orbit coupling. 
Also, the random magnetic field might be parallel to a 
fixed axis in an asymmetric [100] -grown QW due to the 
interplay between the bulk and structural [5(| spin- 
orbit couplings [27! 13 Etjliol l3l| (the structural cou- 
pling depends on the gate voltage [33], so one can tune 
these two couplings to have equal strength). For such 
QW Eq. (J2J is also valid, but in this case unit vector z is 
parallel to the QW plane. In both cases, one component 
of the spin, s z , does not relax. Therefore, these structures 
are especially attractive for spintronics applications. 

In this paper we discuss long-time dynamics of the spin 
polarization in such structures. We consider the relax- 
ation of the vector s = (s x ,s y ) which is perpendicular 
to the random magnetic field (s _L u> p ) and show that 
similar to the velocity autocorrelation function, the spin 
correlation function has long-lived 1/t 2 tail (both for the 
case of strong scatterers and for smooth potential). The 
analogy between velocity and spin relaxation is based on 
the following. As seen from Eq. J2J, the spin rotation an- 
gle is proportional to the integral tp ~ J p x dt and is equal 
to zero for closed paths [13] • Thus when electron returns 
to the impurity its spin restores the original direction. 
This implies some kind of memory effects specific for the 
systems under discussion. 

The Hamiltonian of the system is given by 



p 2 h 



(3) 



where U(r) is a random potential, a z is the Pauli matrix 
and m is the electron effective mass. 

In the Boltzmann approach, the classical dynamics of 
spin related to Hamiltonian Q is described by the kinetic 
equation 



m + (vV)s 



J B S + [Up x s], 



(4) 



where s(r, p) is the spin density related to the averaged 
spin as S = J s(r, p) d 2 rd 2 p/ (2irh) 2 and Jb is the Boltz- 
mann collision integral. Here we consider a case of de- 
generated electron gas (T <C Ep), assuming that the 
spin-polarized electrons have energies close to the Fermi 
energy Ep. First we assume that electrons are scattered 
by strong scatterers randomly distributed in plane with 
average concentration n. In this case 

j B s(9) = nv F [ a(9 - 9')[s(9') - s(6))]d9', (5) 



where a(9) is differential cross-section of one scatterer 
(for electrons with energy E m Ep) and we used short- 
hand notation s(9) = s(r, p) (9 is the angle of the 
vector p). In Eq. J^J we neglected inelastic scatter- 
ing. The role of such scattering will be briefly dis- 
cussed below. To account for classical memory effects 
we will follow the method proposed in Refs. [33] (calcu- 
lation of the velocity correlation function by this method 
is presented in Ref. 0]). The key idea is to replace 
n — > J2i^( r ~ r i) = n + v(r) in the collision integral, 
where z/(r) = ^\ <5(r — rj) — n, (f(r)) = (averaging is 
taken over the position of the impurities). The collision 
integral becomes Jb — ► Jb + J* , where 



J*s(6) = v(r)v J a(9 - 9')[s(9') - s(d)]d9'. (6) 
By the following transformation 



s = f(x)s', f 



cos qx — sin qx 
sin qx cos qx 



q = am (7) 



we eliminate the spin rotation term [u> p x s] from Eq. (0J 



ds' 



+ (vV)s' = (J B + J*)s' 



(8) 



(corresponding unitary transformation of Hamiltonian 
© is presented in Refs. |M El). Following 0, we 
solve equation (J^J treating the term proportional to v(r) 
as a small correction. In the second order of perturbation 
theory we obtain the following equation: 



9s' 



(vV)s' = J B s' + 6Js', 



dt 

5 J = (J*GJ*) 



(9) 
(10) 



where the kernel G(r, tp, ip', t) of the operator G obeys 
dG 



dt 



+ (vV)G = J B G + 5(r)S{p-p')5(t). (11) 



To calculate the average in the Eq. I|10|l we take into 
account that (z^(r)^(r')) = nS(r — r'). As a result we get: 



SJs' = vn / dt'dO' Sa(9-9',t')[s'(9',t-t')-s'(9,t-t% 
Jo 

(12) 

where 



5a(9 -9',t)=v J [a(0 - tp) - a S(6 - tp)} (13) 
x G(0, ip - p', t)[a(<p' - 6') - <7 Q 5(ip' - 0')}dipdp' 

Here cto — J dipa(ip) is the total cross-section and 
G(0, ip — (p',t) — G(r,<p,<p',t)\ T ->Q is the probability for 
an electron starting in the direction n = (cos tp, sin p) to 
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FIG. 1: Four processes giving leading contribution to long 
lived tail in velocity and spin correlation functions. 



return to the initial impurity after time t along the direc- 
tion n' = (cos ip', sin p') (see Fig. QJ. Four terms in the 
product [aid-ip)-a a 6(0-tp)}[(T(tp' -e')-a a 6(tp' -0')} cor- 
respond to four types of correlations [8j shown in Fig. ^ 
Fig.^i shows the process, where electron experiences two 
real scatterings on the same impurity. In Eq. (|13fl the 
corresponding contribution is presented by the term pro- 
portional to a(9 — p)a(p' — 9'). In the process shown in 
Fig. \T]p an electron passes twice the region with the size 
of the order of impurity size without scattering. Since 
the electron "keeps memory" about absence of impu- 
rity at a certain region of space, there exists a corre- 
lation which is accounted for by the term proportional 
to a" 2 5(9 - p)5(p' - 9') in Eq. JT3J). The interpretation 
of the two other terms is based on the fact [8j that in 
the Boltzmann picture, which neglects correlations, the 
following processes are allowed. An electron scatters on 
an impurity and later on passes through the region occu- 
pied by this impurity without a scattering (see Fig. ^J). 
Another process is shown in Fig. ^i. The contributions 
of the terms uq<j(9 — <p)S(ip' — 9') and (Tq5(9 — p)<j(p' — 9') 
in Eq. I|13[l correct the Boltzmann result by substrating 
probabilities of such unphysical events. 

At t 3> t the return probability reads 



G(0,<p-<f/,t) = 



1 



8n 2 Dt 



I 2 

1 nn' 

2Dt 



(14) 



where l^ 1 = nvotx and dtr = f dtpa(tp)(l — cost/?). Inte- 
grating Eq. (|14f> over angles we get 1/AnDt which is the 
probability to return to the initial point with arbitrary 
angle (diffusive return). The second term in Eq. (|14|l is a 
small angle-dependent correction which is responsible for 
the effect under discussion. Indeed, one can see that the 
first term in Eq. 114fl gives zero contribution to Eq. (|13|1 . 



Calculating the contribution of the second term we get 



6<j{6,t) 



of,, cos 9 



(15) 



47r 2 w t 2 

It worth noting that Sa(9, t)dt has a dimension of length 
and can be interpreted as a correction to the scattering 
cross-section due to diffusive returns taking the time ly- 
ing in the interval [t, t + dt]. 

The diffusion equation can be obtained by standard 
means from kinetic equation © with the use of Eq. Q . 
As a result we find the diffusion-like equation for the 
isotropic part of the spin density sq = (,s(p, r, t))e (aver- 
aging is taken over momentum directions): 



dso 
dt 



= -DA s 



na^T f°° s (t - t')dt' 



An 



t' 2 



(16) 



where A = f(x)AT(x)- 1 = (d/dx + qe) 2 + d 2 /dy 2 and 
e is the antisymmetric tensor: es = [e z x s]. Eq. (|16fl de- 
scribes the spin dynamics in the diffusion approximation. 
It simplifies in the homogeneous case: 



ds 
dt 



so 

TS 



Attts 



s (t-t')dt' 

V 2 ' 



(17) 



Here 1/rg = Dq 2 = (ap^) 2 r/2 is the Dyakonov-Perel' 
spin relaxation rate. The initial condition for l|17|) is 
So(0) = Si (we also assume that So(i) = for t < 0). 
Neglecting the second term in the rhs of Eq. i|17[l we 
get the exponential relaxation s (t) = cxp(— i/rs)sj [l5| . 



This solution is valid until exp(— t/rs) 



r TT S /t 2 . For 



larger times, J dt'so(t — t')/t' 2 ss SiTs/t 2 and one can ne- 
glect the term 8sq / dt in Eq. (|17|) . As a result we find that 
the polarization has a long-lived tail 



so(t) 



natr TT S 



"Si 



rSi, 



(18) 



Ant 2 I Ant 2 

which is positive in contrast to Eq. QJ. 

Eq. Ijl8|l was derived for the case of strong scatter- 
ers with low concentration (nd 2 <C 1). The opposite 
limiting case (weak scatterers, nd 2 3> 1) corresponds to 
the smooth random potential with the correlation func- 
tion (U(r)U(r 1 )) = f K q exp[iq(r-r')]d 2 q/(2ir) 2 . In this 
case, the collision integral can be written as a sum of the 
Boltzmann collision integral Jb = (1/t)(<9 2 /dtp 2 ) and 



J* = 



d_ 

dp 



dt[n x f(r)][n x f(r - vi)] - 



(l/2TTm 2 v 3 F ) J °° q 2 K q dq 



rj dip 

(19) 

and f = 



where 1 /r 

-WU(r)/mv F 36]. One can check that (J*) = 0. Sub- 
stituting Eq. into Eq. (ffT77f) . using Eq. (fl4")l and ac- 
counting for two types of correlations |9j we get 



5Js(p) 



2n 2 d 2 



'dt'lG-^t'f^-^ 



dp 2 



-,//,_ ,,d 2 s(p + n,t-t') 



G"(n,t')- 



dp 2 



(20) 
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where G"(<p,t) = d 2 G(0,ip,t)/dip 2 . The calculations 
analogous to the case of strong scattering centers yield 

, . (P TTq 

= 7a IT Si ' ( 21 ) 

In Eqs. (GDI and d = ^2 n 2 q q 3 dq/ n q q 2 dq. 

Above we assumed that S; is homogenous. For slowly 
varying Sj(r), the derived equations relate s (r, t) with 
Sj(r) provided that the spatial scale of inhomogeneity L 
is large compared to ^Dts ~ 1/ma 37]. One can show 
that in the opposite case I -C L <C 1/ma, these equations 
are also valid relating J dr so(r, t) with J dr si(r). 

Let us briefly discuss the role of electron-electron inter- 
action. Such interaction manifests itself both in inelastic 
scattering and in the additional, with respect to the diffu- 
sion process, decay of density fluctuation due to Maxwell 
relaxation. Such a relaxation partially suppresses the 
long-lived tail in velocity correlation function JQ| leading 
to the faster decay [3^: l/t 2 r 2 B /Dt 3 , where re is the 
2D screening length, which coincides with the Bohr ra- 
dius. In contrast to this, the Maxwell relaxation has no 
effect on the spin dynamics because in the classical ap- 
proximation the spin fluctuations are not coupled to the 
charge fluctuations and, as a consequence, do not lead 
to creation of long-range electrical field responsible for 
Maxwell relaxation. 

As for electron-electron collisions, their characteristic 
time T ee is inversely proportional to T 2 . For relatively 
small temperatures, r ee S> r and electron-electron colli- 
sions do not have any effect on the spin relaxation in the 
classical approximation. In the opposite limiting case, 
T ee <C r, electron-electron collisions might suppress spin 
relaxation 39] . The detailed discussion of this case is out 
of the scope of this paper. We believe that l/t 2 depen- 
dence of long-time polarization is also valid for this case, 
while the coefficient in this dependence might change. 

Finally, we compare non-Markovian tail in the 
spin polarization (Eqs. Q18JI. 1(2} ) with the long- 
lived tail induced by weak localization 01 : s o(t) = 
Si(Ts/7rfci?Zt) exp(— t/r<p), where t v is the phase-breaking 
time. At T = 0, when t v = oo, long-time spin dynamics 
is determined by weak localization. However, for T / 
classical memory effects dominate at t 3> t v . 

In conclusion, we developed a theory of long-time spin 
dynamics for 2D system, where spin-orbit-induced mag- 
netic field is parallel to a fixed axis. We showed that in- 
dependently on the type of disorder the non-equilibrium 
spin polarization in such a system decays as l/t 2 (for 
t — ► oo) due to purely classical memory effects. 
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